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Abstract 



d ' We consider m spinless Bosons distributed over I degenerate single- 

H . particle states and interacting through a fc-body random interaction 

with Gaussian probabihty distribution (the Bosonic embedded fe-body 
^ ' ensembles). We address the cases of orthogonal and unitary symme- 

O ■ try in the limit of infinite matrix dimension, attained either as / — > oo 

or as ?n, ^ cxD. We derive an eigenvalue expansion for the second 

moment of the many-body matrix elements of these ensembles. Us- 

^ ' ing properties of this expansion, the supersymmetry technique, and 

the binary correlation method, we show that in the limit / — > oo the 
ensembles have nearly the same spectral properties as the correspond- 
ing Fermionic embedded ensembles. Novel features specific for Bosons 
arise in the dense limit defined as m — > oo with both k and I fixed. 
Here we show that the ensemble is not ergodic, and that the spectral 
fluctuations are not of Wigner-Dyson type. We present numerical re- 
sults for the dense limit using both ensemble unfolding and spectral 
unfolding. These differ strongly, demonstrating the lack of ergodicity 
of the ensemble. Spectral unfolding shows a strong tendency towards 
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picket-fence type spectra. Certain eigenfunctions of individual real- 
izations of the ensemble display Fock-space localization. 
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1 Introduction 

Random-matrix theory (RMT) was introduced by Wigner [|l| to describe 
the complexity of the spectra and wave functions of atomic nuclei. It was 
soon reahzed that RMT could be applied successfully to a large variety of 
systems [@, ^. Examples include atoms, molecules, atomic nuclei and quan- 
tum dots. For such systems, the RMT modeling is not completely realistic 
since all many-body systems are effectively governed by one- and two-body 
forces. This fact motivated work on the two-body random ensembles for 
Fermions (see Refs. Q to ||10[) and inspired the pioneering work of French 
and collaborators on the subject. Especially, it led to the introduction of the 



/c-body embedded ensembles by Mon and French |TT| (see also Ref. [0). In 
the embedded ensembles (EE), stochasticity is generated at the fc-body level. 
The matrix of the fc-body interaction is then taken between ?Ti-particle states 
with k < m. The latter are obtained by distributing m particles over / degen- 
erate single-particle levels. For k < m, the matrix elements of the /c-body 
interaction in the space of m-body states are correlated. The number of inde- 
pendent random variables is much smaller than in RMT. The question arose 
whether these more realistic ensembles yield the same results as RMT. Early 
numerical experiments performed for interacting Fermion systems |^, |^, |^ 
of rather small matrix dimension showed that spectral fluctuation properties 
of the EE agreed with those of RMT. For Bosonic systems, Manfredi | T3[ 



argued that the spectral fluctuations should coincide with the predictions of 
RMT. More recently Patel et al. fl^ performed simulations for a Bosonic 



case with m > I and concluded that the agreement with RMT carries over 
into the so-called dense limit. As far as we know, there are no results on the 
ergodic behavior of the EE for Boson systems. 

Interest in model Hamiltonians with random two-body interactions has 
resurged in recent years in several areas of many-body physics (see Refs. |TB 



to @]), and the question of possible further differences between such mod- 
els and RMT has resurfaced. For instance, Altshuler et al. [O] suggested 



localization in Fock space as a realistic possibility in a system of interacting 
electrons in a random one-body potential. In a different context, calculations 



using random two-body interactions within the nuclear shell-model |21| and 



in the framework of the Interacting Boson Model p3[ led to the observa- 
tion of pairing properties and band structures in the spectra. These "reg- 
ular" features in the low-lying part of the spectrum are robust and seem 
to contradict the expectations based on the canonical ensembles of RMT. 



Recently, three of the present authors found a novel analytical approach to 



the Fermionic embedded ensembles 31, p3l in the limit of infinite matrix 



dimension (/ — > cxd). The main results of this approach are: (i) Within a cer- 
tain parameter range {2k > m), the average spectrum has semicircle shape, 
and the spectral fluctuation properties coincide with those of RMT. (ii) The 
spectral density changes shape at or near 2k = m and becomes Gaussian, 
(iii) In the dilute limit {k <^m <^l) the spectral fluctuations are completely 
uncorrelated (Poissonian). (iv) The spectral fluctuations change gradually 
from Wigner-Dyson to Poissonian. 

In this paper, we extend our analytical treatment to the case of the k~ 
body embedded Gaussian ensembles of random matrices for Bosonic systems. 
Our motivation for doing so is twofold. First, it is instructive to understand 
possible differences arising from symmetry (fully symmetric versus fully an- 
tisymmetric states). To this theoretical argument we add a second and more 
practical consideration: In an atomic trap, Bosonic atoms occupy partly de- 
generate single-particle states. The interaction will lift the degeneracy. A 
random-matrix approach should reveal the generic features of the resulting 
spectrum. 

Our approach is largely patterned after the one for Fermionic systems 
which in turn is based on the invariance properties of the EE. We focus 
attention on the second moment B'^^'> of the random fc-body interaction in 
the Hilbert space of many-body states. We derive an eigenvector expansion 
formally equivalent to the Fermionic one. We exploit this formal equivalence 
and obtain results for the shape of the spectral density and for the spectral 
fluctuations in terms of the parameters fc, m, and / of the problem. All our 
results are obtained in the limit of infinite Hilbert-space dimension. For 
Bosons this limit is realized either by letting / -^ oo (the same limit as for 
Fermions), or by letting m — ^ oo with k and I fixed. This second case is novel 
and has no analogue in the Fermionic case. We refer to this case as to the 
dense limit. The tools at our disposal are the supersymmetry method |3^, |3^ 



and the binary correlation method. Neither of these methods really works in 
the dense limit except to show that standard Wigner-Dyson statistics will 
not apply. Therefore, we resort here to numerical simulations. 

The Bosonic embedded ensembles are defined in such a way that for k = 
m, they coincide with one of the standard ensembles of RMT. Deviations can 
occur only for k < m. In the dilute limit {k <^m <^l) there is no difference 
between Bosons and Fermions. Therefore, we expect that the embedded 
ensembles for Fermions and Bosons are qualitatively similar as long as fc < 



m < I. This is indeed what we find. The physically most interesting and 
novel case is, therefore, the dense limit. Here, we prove analytically that the 
ensemble is not ergodic. Numerical results for the spectral correlations are 
obtained both by ensemble unfolding and by spectral unfolding. In the latter 
case we find a strong tendency of the spectra towards picket-fence behavior. 
We also show that certain eigenf unctions display Fock-space localization. 

The paper is organized as follows. The Bosonic embedded ensemble are 
defined in Section ^ The second moment -B^'^-' is introduced, and an impor- 
tant identity is derived, in Section |^. The eigenvector expansion derived in 
Section ^ is a central analytical result of this work. This expansion is for- 
mally identical to the one found for the Fermionic ensembles. It enables us 
to calculate the low moments of the interaction (Section |) and to apply the 
supersymmetry method (Section^ in a straightforward way. In Section]^ we 
present numerical results. Section ^ contains our conclusions. A preliminary 
account of some of our results was given in Ref. 



2 Definitions and Elementary Facts 

We consider m spinless Bosons in / degenerate single-particle states with 
associated creation and annihilation operators bj and bj where j = 1, . . . , /. 
The Bosons are coupled through a random A;-body interaction Vk{(3), with 
k <m, given by 

^-^ b] . . .b] bj, . . . bj, 

'^'"" ^..£,../" - •' M.:-,.:W(.......) - « 



Kii <i.,<...<ii.<l 



We refer to k as to the rank of the interaction. As in the canonical case, 
we use the labels (3=1 and (5 = 2 for the orthogonal and the unitary 
ensemble, respectively. The matrix element fji,...,jj.;ii,...,ij. of the fc-body in- 
teraction taken between the single-particle states ji, . . . ,jk and ii, . . . ,ik is 
totally symmetric with respect to ji, . . . ,jk and ii, . . . , ik- The elements dif- 
fering in the sequence of indices {ji . . . jk;ii . . . ik} (except for permutations 
of {ji . . . jk} and of {ii . . .ik} and for symmetries specified by (3) are uncorre- 
lated Gaussian-distributed random variables with zero mean and a common 
second moment Vq. Thus, 



The overbar denotes the average over the ensemble. Without loss of general- 
ity we put f Q = 1 in the sequel. The normalization coefficients A/'(ji, . . . ,jk) 
and N'iii, . . . ,ik) in Eq. (|I]) are defined below and are introduced in order 
to obtain for k = m the canonical ensembles of random-matrix theory. The 
number of independent random variables is given by i^/j = P{~^ jT )\.{~^ iT )~^ 
S/3i]/2. This defines the Bosonic fc-body embedded Gaussian orthogonal (uni- 
tary) ensemble of random matrices, respectively, in short BEGOE(fc) and 
BEGUE(fc). 

Hilbert space is spanned by the orthonormal ?7i-particle states |/x) given 
by 

b\ ...h\ 
1/^) = ^ ^ |0) , (3) 

with ji < J2 < • • • < 3m indicating the single-particle levels occupied by 
the m Bosons. Another equivalent definition of the states |/i) is obtained by 
choosing partitions of m, i.e., sequences of positive integers {ni,n2, . . . ,nj.} 
with r < / and X]a=i '^a — '"^; ^^^ ^^^^ ][ < i'2 < . . . < jl of single-particle 
levels. For every such choice, each of the levels j^, s = 1, . . . , r is occupied 
by Us Bosons, and the corresponding ?7i-particle state has the form 

1/^) = \ ,\ |o) . (4) 

The states |/i) are normalized. States belonging to different partitions or, for 
the same partition, to different sequences ][ < i'2 < ■ ■ . < j'^ of single-particle 
levels, are mutually orthogonal. Comparison of Eqs. (§) and @ yields an 
implicit definition of the coefficients A/'(ji, . . . ,jm)- These same coefficients 
appear also in Eq. (0). The dimension of Hilbert space is 

N^C^-'). (5) 

Eq. (g) can be used to define classes of states. All states in a given class 
have in common the same set {ni, . . . ,nj.} of positive integers appearing 
in Eq. (^, irrespective of the choice of single-particle states j[, . . . , j^. By 
definition, permutations of the I single-particle levels leave the class invariant. 
We note that for Fermions there is only one class of states. The number 
A'^ci of classes is given by the number of partitions oi m = ^ n^ subject 
to the constraint rii > n2 > . . . > rir {r < I) . The number of non-zero 



matrix elements coupling the state |/i) to other states depends on the class 
to which l/i) belongs. Therefore, a graphical representation using vertices 
for the Hilbert space vectors and links for the non-zero non-diagonal matrix 
elements (cf. Ref. ||3^) leads to a non-regular graph in general. However, the 
subgraph defined only among states of the same class is a regular graph, i.e., 
the number of links emanating from each vertex is the same for all vertices 
in the subgraph. 

3 The Second Moment of the Interaction 

As in the Fermionic case, the second moment of the matrix (i/|Vfc(/3)|yu) is of 
central importance in our approach. We show that the duality relation found 
for Fermions also holds for the Bosonic ensembles. 

3.1 The Second Moment 

By virtue of the randomness of Vk{P), the elements of the matrix (z/|Vfc(/5)|/i) 
are random variables with a Gaussian probability distribution and zero mean 
value. The distribution is completely specified in terms of the second moment 
Blu\pa{P) defined by 

4t!p. (13) = JiWPMWMF) ■ (6) 

To simplify the notation, we define the operators 

, b\ ...b\ 

and likewise for the adjoint Ba{r)- The index a{r) is a short-hand notation 
for the rank r and for the sequence of indices { ji . . . jV- } • The matrix ele- 
ments {^\B',.Ba(r)W) are real. Using this fact, Eqs. (|l]) and (0), and the 
Hermitecity of Vfc(/3), we obtain 

a{k),-y(k) 



The second moment, Eq. (^, has the same central importance for the theory 



developed here as its counterpart, Eq. (14) of Ref. [^, for the Fermionic en- 
semble. For k = m and (3 = 2, Eq. (|) reduces to 5^u5pa, and correspondingly 
for (3=1. This shows that for k = m, our ensembles reduce to the GUE or 



GOE, respectively, with the radius of the semicircle growing as v iV. 

3.2 Duality 

The duality relation establishes a connection between the second moment of 
the fc-body interaction and that of the {m — A;)-body interaction. The proof 
presented here is somewhat more general than the one given in Ref. [^ . 



We consider Eq. (|^) for the unitary case /? = 2. We observe that the bra- 
state {fi\B',f^s contains (m — k) Bosons and yields the vacuum state upon 

the action of a uniquely defined operator -Bl, _^-,. For the matrix element in 
Eq. (IP not to vanish, the ket-state B^^k)^) niust contain the same (m — k) 
Bosons as {fi\B,f,y In other words, this state yields the vacuum state by 
the action of the operator Bp(^rn-k) with the same index (i{m — k). The 
(m — /i;)-Boson states B\,_^AQi) are orthonormal and complete, i.e., we have 

for fixed (m - k) that E/3(m-fc) '^J(m-fc)|0)(0|'^/3(™-fc) = l(m-fc), with l(.m-k) 
the unit operator for states with {m — k) Bosons. Likewise, the fc-Boson 
states B,^AQi) are orthonormal and complete for the states with k Bosons. 
Hence, 



/3(m-fc),5(m-fc) 



-k) 

-k) 

^wli^-A)\^)m,ik)Bp^m-k)\^) 

5Z (/^l4(m-fe)^<5(--fc)l^)(/^l^l(m-fe)^/3('n-fc)l^) • (9) 



a(fc),7(fc) 
/3(m-fc),(5(m-fc) 



I3(ra-k) 
S{m — k) 



The first equality in Eq. (^) follows from the definition of the second mo- 
ment and the completeness of the {m — fc)-Boson states. In the second one, 

8 



we first rearrange the order of tlie factors conveniently and then use the 
commutation properties of the Bosonic operators. The last line follows from 
the completeness of the fc-Boson states. Using the definition of the second 
moment, Eq. (||), we obtain 

Blll^{2) = Bj::j\2) . (10) 

This is the duality relation. The duality relation is different from and has 
nothing to do with the particle-hole symmetry that is used for treating 
Fermions in more than half-filled shells. 



4 Eigenvector Expansion of the Second 
Moment 

In the Fermionic case, results for the shape of the spectrum and spectral 
fluctuations could be derived with the help of the eigenvector decomposition 
of the second moment. A similar decomposition exists in the case of Bosons 
and is derived below. A central role in the derivation is played by the identity 

V- bj,...bj^bl...bl ]_ Y^ t ^t ni) 

This identity relates a sum over the indices ji, . . . , jfc restricted by the con- 
dition ji < . . . < ik containing the appropriate weights to the unrestricted 
sum over the same indices. The weight factor fc!/[A/'(ji, . . . , Ja:)] is related 
to the number of permutations with repetitions of k indices. A proof of this 
identity is given in the Appendix. 

We consider flrst the unitary ensemble. Proceeding as in the case of 

(fc) ^{sa) _ \{s)f7\Msa) 



Fermions, we solve the eigenvalue equation ^p^B^J^paCap = A^'{k)Ci 



flU 



and show that the eigenvectors have the form 

Clf = {a\Bl^^^Bsis)\p) . (12) 

Here, s is the rank of -BL n and Bs^s), and a labels the pair {l3{s),5{s)). We 
assume first that no index in f3{s) equals any of the indices in S{s). Using 
Eqs. (P) and (|T2|) in the eigenvalue equation, we obtain 



(.) (2)c(-)= V (/i|SL,,S,(.)S,(,)St 4, .S„(,)|^) . (13) 



(XT a{k)'y{k) 



We have used that there are no common indices in (3{s) and S{s), so -BL ^ and 
Bs{s) commute. To perform the summation over 7(A;) we use the identity (|lTl) 

and the fact that ^- bjb] = l + Y.j b]bj. This yields (^+™r~^)- ^or the sum 
over a{k) we proceed similarly, rearranging first the order of the operators 
to bring together the operators involving a{k). We use an identity similar 
to ( p!TD to compute the summation. We obtain the factor (™^'*) multiplying 

Cjfu ■ This shows that the matrix Cap defined in Eq. (|12]) is indeed an 
eigenvector with eigenvalue 

(s) / , N /'^ — s\fl + 171 + 3— 1 



In Eq. (|Tip and below, the label B stands for Bosons; we shall use the label F 
for Fermions. The eigenvalues Ag (A;) decrease monotonically with increasing 
s and terminate at s = m — k, as for the case of Fermions. Again, this is 
consistent with the GUE result for k = m, where only the term s = occurs 
and where Ag (m) = N. 

Eigenvectors of the form (^) with no common index in j3{s) and 6{s) 
are orthogonal. In the case where there is a common index, a construction 
identical to the Fermionic case is required. The eigenvectors can be "lifted" 
from a lower rank s' to a higher one with s > s' by insertion of powers of 
the number operator. The construction of linearly independent eigenvectors 
(j(sa) proceeds as for Fermions (see Ref. [Q for details), and one finds that 



these eigenvectors posses the same eigenvalues (|14|). As in the Fermionic 
case, the eigenvectors do not depend upon the rank k of the interaction, 
only the eigenvalues do. The construction of the eigenvectors is identical 
in both cases. The dimension D^ of the subspace spanned by degenerate 
eigenvectors characterized by s is given by D^ = 1 and for s > 1 by 

Js) fl + s-lV fl + s-2 " 



It follows that Xir^o-^B ~ ^^- Therefore, the eigenvectors form a com- 
plete basis. We impose Hermitecity on the eigenvectors and normalize them 
according to 

Y.'^i'^^^^^f = ^^^^^^'^ ■ (16) 

flU 

It is instructive to compare the eigenvalues for Fermions and for Bosons. 



We recall that for Fermions, the eigenvalues have a form similar to Eq. (14) 

10 



with the second factor on the right-hand side replaced by ( ^™^ ^^) . Since 
the dimensions of the Hilbert spaces for Bosons and for Fermions differ, and 
since the eigenvalues depend on the dimension, a fair comparison is obtained 
by normalizing the eigenvalues to the value A^^^^ (k) . For Bosons and Fermions, 
the ratios i^ (k) = A^ {k)/Al (k) with i = F, B are given by 



4^) (A:) 



?^^\k) 



n 

k 

n 



m + 1 — i — s 
m + 1 — i 

m + 1 — i — s 
m + 1 — i 



I -\- m — i + s 
I + m — i 

I — m + i — s' 
I — m + i 



(17) 
(18) 



(«)/ 



rM)l 



^(0), 



It is easy to see from Eqs. (0) and (|T|) that tp'{k) < e^'{k) while tp'{k) = 
£g (/c) = 1 by definition. In the dilute limit k <^m <^l where the distinction 
between Fermions and Bosons is irrelevant, one finds the expected result 
£p (fc) ~ £g (A;). On the other hand, for fixed k and /, the difference between 
the Fermionic and the Bosonic case increases with increasing particle number. 
For Bosons, the eigenvalues with larger s become more important. This fact 
suggests stronger deviations from Wigner-Dyson statistics for Bosons than 
for Fermions. This is particularly obvious for values of m beyond half-filling. 



W/ 



9^/ 



2m > /, and m — k > I — m, where ip (k) = ioi s > I — m while ^b (A;) 
is different from zero. This is demonstrated in Fig. |l] which shows i^^\k) for 
Bosons and for Fermions as a function of s, for k = 2, I = 40 and for several 
values of m. 

For the orthogonal case we observe that the matrix B^''\l) in Eq. (H) 
consists of two terms, each one of the form of B^''\2). Combining the ex- 
pressions for (3=1 and for /3 = 2, we find for the second moment the 
eigenvector expansion 



5(^) 



fiu^pa \ 



(^) = l^T.^Bik) 



^ in, ^ nrr I '^ P 1-^ fj,p ^ U(J 



flU 



pa 



(19) 



This expression shows a formal analogy to the Fermionic case. However, the 
dimension of Hilbert space, the values of the Ag (/c), and the definitions of 
the eigenvectors, are all different. 
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Figure 1: The ratio i^^\k) = A^'^\k) / A^^\k) as a function of s for (a) Bosons 
and (b) Fermions, for k = 2, I = 40, and for several values of m. Only 
non-zero eigenvalues are shown. We note the difference in horizontal scales 
in the two plots. 



5 Low Moments of Vk 

We apply the results obtained in Sections ^ and ^ to calculate the dependence 
on k, m, and / of measures of the shape of the average spectrum in the 
limit N ^ oo. Only in this limit can we expect to obtain generic results. 
In contradistinction to the Fermionic case, where the limit N —>■ oo implies 
/ ^ oo, the limit N -^ oo can for Bosons also be obtained by letting m ^ oo. 
In particular, we will be interested in this case for fixed values of / and k. We 
refer to this case as to the dense limit. Having chosen the Bosonic ensembles 
in such a way that they agree with GUE or GOE for k = m, and observing 
that for k <^m <^l Bosonic and Fermionic ensembles must agree, we do not 
expect a qualitatively different behavior for Bosonic and Fermionic systems 
as long as TTi < /. Thus, the dense limit is the really interesting and novel 
case for Bosons. It has no analogue for Fermions. 

We consider three ratios that yield information about the shape of the 
average spectrum. The ratio S measures the fluctuations of the center of the 
spectrum in units of the average width of the spectrum. The ratio R measures 
the relative fluctuation of the width of the spectrum. These two quantities 
yield information about finite-size effects in numerical calculations [Q . The 
third ratio is the kurtosis k, written as k = Q -|- 2. The quantity Q marks 
the difference between the semicircular {Q = 0) and the Gaussian shape 
{Q = 1) of the average spectrum. The definitions of S, R, and Q are given in 
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Ref. |32|. Since the eigenvectors and the duahty relation are formally similar 
to the Fermionic case, the formal expressions for the three quantities S, R, 
and Q are identical, too. Differences are due to the changed form of the 
eigenvalues, to different ranges of the parameters k, m, and /, and to the 
following fact. 

The creation and annihilation operators for Bosons are not nilpotent like 
those for Fermions. Therefore, contributions of the form '^ C^v Cj^p that 
for Fermions are of relative order (9(1//^) and vanish in the large N limit, 
may become important for the Bosonic ensembles. Such terms can be worked 
out explicitly exploiting the symmetry or antisymmetry of the Hermitean 
eigenvectors C\^v of the second moment. We define the "parity" TXa{s) as 
zero (one) when C\^v is symmetric (antisymmetric) in the indices /x, i/, and 
have C/j'T = (— 1 )'""'"' Ci/^". Using duality and the orthogonality relations, 
Eqs. (|10|) and 

E.„(-i)^"'=' = 

important identity for the partial sum rfg'' = X]a(~-'^)'^" 



m 



^^^A('')(m) ^ C/tt C^tt , we obtain the sum rule 



N . From this sum rule, it is easy to derive the following 



d 



(s) 



B 



1 + S 



l + s-2 
s-1 



l + s-2 
s 



(20) 



Using Eq. (| 
S{k,m, I) 



and the result for the Fermionic ensembles, we obtain 

-1 



,(0) 



l + 6^MBi^-k) tT[Vk{f3)] 



(21) 



r '^ 

R{k,m,l) = 2Ul + 6^,)J2i^^BfDB+^hiY.(^B\m-k)fd^^' 

'^ s s=0 



X 



tr[Vkm 



(22) 



and 



, min(m— fe,A;) 



Q{k,mJ) = ^1 E ^^B\k)A^B\^-mD^B^ + ^hid^B^] 

m—k „ „ 

+ V [E «"(*))"< -f(I^ ''="<) 

s=0 s 
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{sa),{tb) 
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X 



j^tr[V, 



(23) 



The width of the spectrum is given by 



1 

N 



trWMV=AS\k) + '-^j:A 



B "B 



(24) 



In Eqs. (pi|)-(P^), traces are always taken in the Hilbert space of the many- 
body states \fi). We have simphfied the notation by omitting the dependence 
on A; or m — A; of the eigenvalues Ag in the terms which are invariant under 
the exchange k ^^ m — k because of the duality relation. 

We consider the large A^ limit, first realized by letting I ^ oo, with m/l 
fixed. The results are similar to the Fermionic case: For k > 1 both S ^ 
and R ^ 0. The convergence to the limit is particularly slow when both k 
and the filling factor m/l are fixed. In this case, S oc l~^ oc (IniV)"'^ and 
R oc l^"^^ oc (InA^)"^'^ while for the canonical ensembles {k = m) the limit is 
approached as S{m, m, I) oc 1/A^^ and R{m, m, I) oc l/N"^. The behavior of 
Q suggests that the shape of the spectrum coincides with the semicircle for 
2/c > m. The semicircle undergoes a smooth transition to Gaussian shape for 
2k < m. The critical value of k where this transition takes place is associated, 
as in the case for Fermions, to duality. 

We turn to the dense limit, letting m -^ oo with both k and / fixed. We 
find 



lim S{k,m,l) 



2k\ fl+k-l\-^ 



(i+^.i)(.)rr^) 



f2k 



I') + ^/^i E 



( 2k \ n+k+s- 



s=0 \k+s 



k+s 



y'd 



is) 



(25) 



lim R{k,m,l) 






(26) 



lim Q{k, 171,1) 



1 



2k 
k 



-1 



2k 

k + s 



l + k + s- 
k + s 



-1 



D 



(s) 



(27) 
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The identity used to obtain the second equahty in Eq. ( P7[ ) is derived 
in Ref. [0. Eq. (^) apphes only in the unitary case (/? = 2) and imphes 
that in the dense hmit, the average spectrum has Gaussian shape. We have 
been unable to extend this result to the orthogonal case (/3 = 1). This is 
because we did not succeed in calculating the term tr(C*^'^"''C*^*"''C''*''^C^*^^) 
appearing in Eq. (p3|). On physical grounds, however, we expect an equation 
analogous to Eq. (|27| ) to be valid also for the BEGOE(A;). We conclude that 
the shape of the average spectrum has Gaussian form. This is in keeping 



with the results of Ref. 38 



A more important and surprising result lies in the fact that the right-hand 
sides of Eqs. (|^) and (|^) do not vanish: The fluctuations of the centroids 
and of the variances of individual spectra do not vanish asymptotically. This 
feature differs from the behavior both of canonical random-matrix theory 
and of the embedded Fermionic ensembles. We are led to the important 
conclusion that the Bosonic ensembles are not ergodic in the dense limit m — »• 
cxD with k and / fixed: Unfolding the spectra by taking the ensemble average, 
or the spectral average over a single member of the ensemble, will yield 
different results. We attribute this non-ergodic behavior of the ensemble 
to the slow rate at which the eigenvalues Ag (fc) decrease as s increases. 
Physically, it is a consequence of the fact that the number of independent 
random variables Kp in the ensemble does not grow with m but stays fixed 
with k and /. We return to this point in Section |^. 

6 Supersymmetry Approach 

As in the case of Fermions, the eigenvector expansion (|19]) permits us to 
use the supersymmetry approach. We shall not give a self-contained pre- 
sentation of this technique in the present context. There are review papers 
which describe it extensively. We mention for instance, Refs. |^^ and pl |. 



Moreover, the supersymmetry technique was also described in some detail 



in our preceding paper on Fermions |52[. Suffice it to say that the result- 



ing expressions are formally the same as in the Fermionic case. Because of 
specific properties of the eigenvalues Ag (fc), in particular in the dense limit, 
new contributions appear in the lowest-order terms in the loop expansion. 
These contributions may be finite and, thus, imply departures from RMT 
behavior. In the dense limit, such behavior must be expected because of the 
non-ergodic features found in Section ^ 
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6.1 Saddle Point 

We summarize the main steps of the supersymmetric approach for the one- 
and two-point functions; details may be found in Ref. [^. After averag- 
ing over the ensemble, the integrand of the generating functional contains 
an exponential whose argument is a sum of squares of bilinear forms in the 
integration variables involving the matrix B^''\ We use the eigenvector ex- 
pansion ( p!9D and perform the Hubbard-Stratonovich transformation. For 
each eigenvector C^^'^^ , this introduces a supermatrix a'-'^"-' of composite vari- 
ables. We use the saddle-point approximation. We define E = {Ei + E2)/2 
and e = E2 — El, where Ei and E2 are the energy arguments of the advanced 
and retarded Green functions, respectively. Omitting small terms of order 
1/A^, we find the coupled saddle-point equations 



^isa) 



^tr,([E - 5^ Ai*)(A:)a(*^)c(*^)]-^A«(A:)C(-)) , (28) 



N 

tb 



where Ag (fc) is the positive square root of A^ (k). As in the Fermionic 
case, this set of equations is solved iteratively by introducing the variable 
Xpa = Esa Ai3'^(fc)a(^'^)cir^ and using the identity G{E) = [E - X]-^ for the 
averaged Green function. We find that X obeys the saddle-point equation 
X = A^ (k)[E — X]~^. Hence, the solution diagonal in the superindices has 

the form X;|jfg = (5^^AS5'^(fc)r(°), with 



for l-E] < 2Ag . The ± signs refers to the retarded {G~{E)) and advanced 
{G~^{E)) case, respectively. This implies that the solutions of Eq. (pSD reduce 
to r*^°) for s = 0, while a^^""^ = for s > 1. For the one-point function this 
shows that, within the range of validity of the saddle-point approximation, 
the semicircle describes the shape of the average spectrum. 

As in the case of the canonical ensembles, the invariance of the effective 
Lagrangean under general pseudounitary transformations implies that the 
two-point function possesses not a single saddle point but a saddle-point 
manifold which is generated by T^^a^^^T. The matrices T parameterize the 
coset space U0SP(1, 1/1, 1)/[U0SP(1, 1)®U0SP(1, 1)], and a(°) is a diagonal 



supermatrix of dimension four with entries given by Eq. (2£) in the usual 
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way EM. Using this and the saddle-point solution in the expression of the 



effective Lagrangean, the first-order term in e takes the canonical form 

77rf 

- — trg(LT-iLT) . (30) 

Here d is the average level spacing and the diagonal supermatrix L (defined in 
Ref. ||3^) distinguishes the retarded and advanced cases. The argument also 
carries through for higher-point correlation functions. Eq. (^0|) shows that 
within the range of validity of the saddle-point approximation, the spectral 
fluctuations of the BEGUE(fc) are identical to those of the GUE. 

6.2 Corrections of the Loop Expansion 

To determine the range of validity of the saddle-point solution in the limit 
N ^ oo, we use the loop expansion, a power-series expansion of the effec- 
tive Lagrangean around the saddle-point, followed by an expansion of the 
resulting exponential. Again, the formal equivalence between the Fermionic 
and Bosonic eigenvector expansions permits a straightforward transcription 
of the results for Fermions. 

For the one-point function, the first correction term in the loop expansion 
yields 

\Gm\<^^^^U + Q{k,rn,l)) . (31) 

From Eq. ( ^1]) we conclude that for I —>■ oo, the shape of the spectrum is the 
semicircle for 2k > m, and finite corrections arise for 2k < m (see Section ^. 
These results coincide with those for the Fermionic ensembles. In particular, 
in the dilute limit a Gaussian shape of the spectrum is obtained. This holds 
also true in the dense limit m -^ oo with k and / fixed as given by Eq. (pT]). 
For the two-point function we drop terms that vanish as 1/A^ for A^ — > oo. 
The resulting expression can be written as an exponential which contains two 
non-trivial terms. The first one contains the expression 

- ^trg(LT-iLT) , (32) 

"1 

where di is a modified level spacing which includes the terms that lead to the 
bound given in Eq. (|3l| ) for the average Green function. More interesting is 
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the second term which, as in the Fermionic case, is proportional to 

R{k,m,l)(^tTg[-eLTLT-^ + LTJT-^]) . (33) 



Kravtsov and Mirhn [^ have shown that terms of this form produce devia- 
tions from universal behavior. The strength of the fluctuations which arise 
from this term is proportional to the ratio R{k,m,l) given in Eq. (0). We 
recall that R measures the relative fluctuations of the width of the spectrum. 
While in the dilute limit this quantity tends slowly to zero as iV — > cxd, in 
the dense limit it attains a constant value. Such finite corrections imply that 
the saddle-point solution is no longer appropriate in the dense limit, and 
deviations from Wigner-Dyson statistics must occur. 

We conclude that in the dense limit, BEGUE(A;) possesses spectral fluctu- 
ations which deviate significantly from the predictions of RMT. As mentioned 
above, we attribute this feature to the fact that the number of independent 
random variables of the ensemble remains the same as m increases, see Sec- 
tion 1^. Alternatively, we attribute this result to the rather slow decrease of 
the eigenvalues Ag with increasing s. Contributions other than the canon- 
ical term s = are relevant and give rise to important corrections for the 
spectral fluctuations. These features are similar to those characterizing the 
dilute limit. Here, the spectral correlations become Poissonian. In that case, 
however, the function R{k, m, I) vanishes asymptotically, albeit very slowly. 
Therefore, the first loop correction of the supersymmetry method is not able 
to provide solid evidence for non-RMT spectral fluctuations in the dilute 
limit. However, in this limit another method is available: The extension of 



the binary correlation approximation of Mon and French |]TT| developed in 



Ref. |3^. This method works equally for Bosons and for Fermions (in the 
dilute limit, the distinction between the two cases disappears) and, thus, 
yields identical results: The spectral density acquires Gaussian shape, and 
the spectral fluctuations become Poissonian. 

7 Numerical Results: the Dense Limit 

In the dense limit, there are no analytical techniques to predict spectral 
fluctuations. Moreover, the Bosonic ensembles are not ergodic. Therefore, 
unfolding procedures based upon averaging over the ensemble carry a ques- 
tion mark and spectral unfolding is the preferred procedure to study the 



spectral fluctuations. In this Section, we present numerical results on the 
spectral fluctuation properties of BEGOE(fc) for k = 2 in the dense limit. 
We exhibit the difference between spectral unfolding and unfolding by us- 
ing the average over the ensemble. For brevity, the latter method is called 
ensemble unfolding. 

Numerically, the dense limit is conveniently modeled by taking 1 = 2. The 
dimension of Hilbert space is then given by A^ = r?7, + 1 (cf. Eq. (^), and the 
situation m ^ / can be easily attained numerically. Furthermore, for / = 2 
and k = 2 the number of independent random variables of the ensemble is 
Ki = 6. We emphasize that Ki is independent of the number m of Bosons. 
The many-particle states |/i„) are written as {m — n,n) where m — n and 
n indicate the number of Bosons occupying the flrst and the second single- 
particle state, respectively. The many-particle states |/i„) are arranged in a 
sequence (m, 0), (m — 1, 1), . . . , (0, m) of increasing values of n. Then, the 
Hamiltonian matrix attains a band structure, with non-zero matrix elements 
on the main diagonal and on the k closest diagonals. We restrict ourselves 
to the physically most relevant case k = 2. 

7.1 Ensemble Unfolding versus Spectral Unfolding 

The evaluation of measures of spectral fluctuations requires that individual 
spectra be unfolded: The average level spacing must be constant. Spectra 
can be unfolded in two ways. Spectral unfolding consists in handling each 
spectrum independently. One looks for a transformation which yields the 
value unity for the average level spacing. This transformation need not be 
the same for different spectra. Ensemble unfolding uses the average spectral 
density (obtained by averaging all spectra in the ensemble) to flnd a single 
transformation which is then applied to all spectra of the ensemble. 

Spectral unfolding is typically done by fltting a polynomial to the staircase 
function of each realization of the ensemble. In the results presented below, 
we have determined the degree of the polynomial for each realization by 
taking the value which gives the flrst minimum of the associated x^ when 
increasing the degree starting from 1. The maximum degree considered was 
20. Typically the flrst minimum was found around 11. 

Ensemble unfolding is usually preceded by removing the fluctuations of 



the center and of the width [T^, |36|. All spectra are recentered and rescaled 
so as to have a common center and width. This procedure is essentially 
equivalent to the more sophisticated method used in Refs. 0, 0, which 
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Figure 2: Spectral density of the BEG0E(2) for / = 2 and m = 3000, 
normalized to the dimension A^ = 3001 of Hilbert space. The curves (a) and 
(b) show the level densities of two individual realizations of the ensemble; 
curve (c) shows the ensemble-averaged spectral density. The dashed curve 
corresponds to a Gaussian shape of the spectrum. 

corrects higher moments of the distribution and relies strongly on the Gaus- 
sian shape of the average spectrum. The recentered and rescaled individual 
spectra are then superposed to obtain the ensemble-averaged spectral den- 
sity and the staircase function for the ensemble. The latter is fitted with a 
polynomial, and each spectrum of the ensemble is unfolded using this polyno- 
mial. However, for the Bosonic ensembles in the dense limit, this procedure 
yields a non-Gaussian average level density as a consequence of the non- 
ergodic behavior of the ensemble. Therefore, we did not recenter and rescale 
the individual spectra for constructing the ensemble-averaged spectral den- 
sity and the corresponding staircase function. We have used a polynomial 
of degree 11 for fitting the average staircase function. The stability of the 
results has been checked by varying the degree of the polynomial and the 
size of the bins used to construct the average spectral density. 

Fig. ^ shows the spectral densities of two realizations of the ensemble 
and the ensemble-averaged spectral density obtained from 1512 realizations 
of the ensemble for m = 3000. The integral over each histogram is normal- 
ized to 3001, the dimension of Hilbert space. In the Figure, we have also 
plotted a Gaussian density (dashed curve), whose width is given by the the- 



20 



oretical prediction Eq. (P^). Comparison of the ensemble-averaged density 
with this curve confirms the prediction of a Gaussian shape of the average 
spectrum. The striking differences between each of the two spectral densities 
and between those and the ensemble-averaged spectral density illustrate the 
non-ergodic character of the Bosonic ensemble for m ^ I. The Gaussian 
form of the average spectrum arises as a consequence of the Central Limit 
Theorem. It has no bearing on individual spectral shapes. 

7.2 Spectral Statistics 

We turn to the spectral fluctuation properties of BEG0E(2). We have ana- 
lyzed 1512 members of the ensemble for m = 3000. After fitting the polyno- 
mials, including all levels of each realization for both spectral and ensemble 
unfolding, we only considered 60% of all levels, i.e., 1800 levels for each 
realization. These were the levels closest to the center of the spectrum. 

In Fig. ^ we present the nearest-neighbor spacing distribution P{s) ob- 
tained by ensemble unfolding and by spectral unfolding. In both cases the 
distribution P{s) corresponds neither to the Wigner surmise nor to a Poisson 
distribution. In the case of ensemble unfolding the level repulsion character- 
istic of the GOE is clearly lost. On the other hand, the tail of the distribution 
decays rather rapidly in comparison to the exponential decay of the Poisson 
distribution. For spectral unfolding the distribution P{s) is dominated by 
a prominent peak centered at s = 1. This suggests that individual spectra 
have an almost constant level spacing, i.e., are nearly of picket -fence type 
(one-dimensional harmonic oscillator). 

The two-point correlations of the spectra are characterized by the A3- 
statistics of Dyson and Mehta measured at the centers of the spectra and 
shown in Fig. ^. We have computed A^{L) directly from the spectra using 
the expression derived by Bohigas and Giannoni [^. For comparison we have 
also plotted the results for a GOE spectrum A3(L) ~ [ln(27rL) + 7 — tt^/S — 
5/4]/7r^, an uncorrelated (Poisson) spectrum A3(L) = L/15, and a picket- 



fence spectrum A3(L) = 1/12 (see Ref. |^ for details). We have checked that 



the increase of A3(L) with L for small values of L is approximately linear in 
both panels of Fig. ^ Ensemble unfolding (Fig. ^(a)) yields a A3-statistics 
which deviates rapidly from GOE behavior, and increases almost linearly. 
This is consistent with the expected deviations from RMT predictions found 
in Section |[ In the case of spectral unfolding, A3(L) is almost constant up 
to L ~ 20. This result again suggests a tendency of individual spectra to 
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Figure 3: Nearest-neighbor spacing distribution P{s) obtained by (a) en- 
semble unfolding and (b) spectral unfolding. The dotted curves correspond 
to the Wigner surmise and the dashed ones to the Poisson distribution. We 
note the different vertical scales used in the frames. 

have a picket-fence-like behavior. Beyond this point, the A3(L) grows albeit 
less rapidly than after ensemble unfolding. 

7.3 Structure of Individual Spectra and Wave Func- 
tions 

The results presented above show the deviations from RMT behavior ex- 
pected after ensemble unfolding. For spectral unfolding the results suggest 
that spectra of individual realizations have a very rigid structure with an 
almost constant spacing in some regions of the spectrum. We now describe 
in more detail the structure of individual spectra. While we do not have a 
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Figure 4: The Aa-statistics (solid lines) measured at the centers of the spec- 
tra after (a) ensemble unfolding and (b) spectral unfolding. For comparison 
we have plotted the results for the GOE (dotted curve), for a Poissonian spec- 
trum (dashed line), and for a picket-fence spectrum (dotted-dashed line). 

full analytical understanding for the tendency towards picket-fence behavior, 
we believe that the following observations are related to the small number 
of independent random variables of the Bosonic ensemble for m ^ I, and to 
the associated graph structure. 

Fig. ^ shows the staircase function 1]{E) for a typical realization of the 
ensemble {k = 2, I = 2) for m = 1000. The spectrum is dominated by levels 
with almost constant spacing. The staircase function typically displays one 
or more points where an abrupt change in the density of states takes place: 
The spectrum has (almost) constant spacing up to the level 619. This behav- 
ior suddenly changes at the level 620. While the spacing between neighboring 
levels is no longer constant, the spacing between next-to-nearest-neighbors 
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Figure 5: Detail of the staircase function ri{E) for a typical realization of 
BEG0E(2) with m = 1000 and I = 2. The vertical lines help to show how 
two (almost) equidistant spectra overlap around the level 620. We note the 
change in the density of states after this level. 

is and is almost the same as it was before level 620 appeared. This is illus- 
trated by the vertical lines plotted above and below the staircase function. 
We remark that these observations imply non-stationary properties of the 
spectra. 

These results suggest that the spectrum of an individual realization con- 
sists of pieces of overlapping segments of spectra with almost constant level 
spacings. The kink at level 620 in Fig. ^ marks the left edge of the overlap 
region. Such spectra would lead to results like those in Figs. |^(b) and ||(b). 
We notice that if we were to use for spectral unfolding a piecewise polyno- 
mial fit instead of a polynomial fit covering the whole spectrum, the rigid 
features of the spectrum would appear in an even more convincing form in 
the spectral statistics. We have not pursued this point further. 

If the spectra do possess the structure proposed above, eigenf unctions 
belonging to different segments of the spectrum should differ qualitatively. 
To verify this hypothesis, we write the eigenvectors \i) as linear combinations 
of the ordered many-body basis states |/i„). 



P/ / J CinlUn) 



(34) 
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Figure 6: Probabilities |cj„p for eigenvectors belonging to eigenvalues in 
the vicinity of the kink shown in Fig. ^. The two overlapping segments 
of nearly equidistant levels are easily distinguished by the structure of the 
eigenfunctions. 
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In Fig. ^we plot the probabilities |cj„p of eigenvectors \i) belonging to eigen- 
values in the vicinity of the kink shown in Fig. ^. Eigenvectors up to i = 619 
behave similarly and are somewhat extended, although the distribution of 
the intensities clearly deviates from the Porter-Thomas distribution expected 
from RMT. This behavior of the eigenfunctions changes abruptly at i = 620 
where the eigenvector is a rather "coherent" superposition of close-lying ba- 
sis states. We emphasize the difference in the horizontal scales used in Fig. ^ 
for i = 619 and i = 620, for instance. The eigenvectors with i = 621,623 
display the same behavior as the eigenvectors up to i = 619, and thus corre- 
spond to the first segment of the equidistant spectra. The eigenvectors with 
i = 620, 622, 624 are much more localized in Fock space. They differ in the 
number of intensity oscillations. As i increases, so does this number, and the 
levels on this sequence become more and more delocalized. At some point 
the spread of the eigenfunctions in the second segment is indistinguishable 
from that in the first one. 

8 Conclusions 

We have studied the shape of the average spectrum and the spectral fluctua- 
tions of the Bosonic embedded ensembles BEGOE(A;) and BEGUE(fc) in the 
limit of infinite matrix dimension, attained by letting either the number / of 
degenerate single-particle states go to infinity, or by letting the number m of 
Bosons go to infinity. In the first case, the results are qualitatively the same 
as for the Fermionic embedded ensembles [^ : For sufficiently high rank k of 



the random interaction {2k > m), these ensembles behave generically. The 
average spectrum has semicircle shape, and the spectral fluctuations obey 
Wigner-Dyson statistics. A smooth transition to a different regime takes 
place at or near 2k = m. This new regime is characterized by a Gaus- 
sian spectral shape. The spectral fluctuations also change and are not of 
Wigner-Dyson type for 2k < m; in the dilute limit k <^ m <^l the spectral 
fluctuations become Poissonian. We expect that the transition away from 
Wigner-Dyson statistics is more rapid in the Bosonic than in the Fermionic 
case. 

In the second case {m > I), and in particular in the dense limit, where / 
and the rank of the interaction k remain flxed asm ^ oo, we have shown that 
the Bosonic ensembles are not ergodic. To the best of our knowledge, this is 
the flrst case that a random-matrix model has been shown to be non-ergodic 
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in the limit of infinite matrix dimension. Moreover, we have shown that the 
spectral fluctuations deviate strongly from RMT results. This finding dis- 



agrees with the conclusions based on numerical simulations of Refs. |T^, |n 
We ascribe this to the fact that the ratio m/l for the parameter settings 
investigated both in Ref. [[ll (m = 11, / = 4) and Ref. (m = 10, / = 5) 



were still too small to encounter deviations from Wigner-Dyson statistics. 
We have studied numerically the structure of the spectra of individual real- 
izations of the ensemble in the dense limit using the case m = 3000, / = 2 
and A; = 2 as example. In both short- and long-range correlations we found 
a tendency towards a picket-fence type of behavior. We presented evidence 
pointing to a non-stationary behavior of the spectra. The spectrum of an 
individual realization seems generically to consist of segments of spectra of 
picket -fence type. These segments overlap. In the overlap region, the eigen- 
vectors can be attributed to the different segments by their structure which 
is characteristic for each segment. In particular, we have shown that some 
of the eigenfunctions display strong localization in Fock space. We cannot 
offer an analytical explanation for these properties. We believe that they are 
caused by the small and constant (independent of m) number of independent 
random variables Kf3 and by the specific graph structure associated with the 
Hamiltonian matrix of the Bosonic ensembles. This causes the Hamiltonian 
matrix to attain a particular structure with highly correlated matrix ele- 
ments. In particular, we do not know how the spectra change as / increases 
for m/l ^ 1. 
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Appendix: Proof of the Identity Eq. ( [tt] ) 

We claim that 



is equal to 



h- h- h^ h^ 
r, = l Yl bn---b,A---bl.- (36) 

jl,J2---,jk 
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We prove the claim by induction. The claim is trivial for A; = 1 and k = 2. 
We assume that the assertion holds for k and prove it for k + 1. On the 
right-hand side of Eq. (pSf), written for k + 1, we use the definition of the 
coefficients J\f{ji, ■ ■ ■ ,jk+i)- We write 

^'^' = ^ ^jTT^g ' ^^^^ 

where the sum over the non-zero integers n, is restricted by Xli^j = k + 1. 
Moreover, j[ < ... < j'^.. We arrange the terms in the sum in Eq. (^) in 
increasing powers of bji (and b-, ) and apply the induction hypothesis to the 

remaining sums involving at most k factors b-. Hence, 
fc+i . 

'^^+1 = Z2 fk + l-n)\n\ 

n=l ^ ' 

Tib] Tbl ...b\ 



^ ik + l-n)\n\ 

n=l ^ ' 




X Y, b,.: 


■ bjk + l-.n\'^jk + 2 


jl,---dk + l-n<jk + 2-n 




1 y^A + l^ 




X Y ^h ■ ■ 


■ "ifc+l_„l"ife + 2 


jl,---dk + l~-n<jk+2-n 





"(b] Vbl ...b\ . 

^ Jk+2-n' 3k+l-n Jl 

(38) 

Consider now the sum on the right-hand side of Eq. (^) for k + 1. For 
each term, we identify the largest of the indices ji, . . -jk+i- We call it jmax- 
The ftj^ax's and bj^^Js come in powers rij^^^^ = 1, . . . , fc + 1. For a fixed value 
^imax = ''^5 the number of terms in the sum is ( ^ ) • Relabeling the indices 
in the sum, so that jmax -^ jk+2-n, the sum can be rewritten in the form of 
Eq. (H). 
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